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CP^ 1 model with the theta term and maximum entropy method * ' * 

Masahiro Imachi a , Yasuhiko Shinno b and Hiroshi Yoneyama c 
a Department of Physics, Yamagata University, Yamagata 990-8560, Japan 
b Graduate School of Science and Engineering, Saga University, Saga 840-8502, Japan 
c Dcpartment of Physics, Saga University, Saga 840-8502, Japan 

A 8 term in lattice field theory causes the sign problem in Monte Carlo simulations. This problem can be 
circumvented by Fourier-transforming the topological charge distribution P(Q). This strategy, however, has a 
limitation, because errors of P(Q) prevent one from calculating the partition function Z(8) properly for large 
volumes. This is called flattening. As an alternative approach to the Fourier method, we utilize the maximum 
entropy method (MEM) to calculate Z(8). We apply the MEM to Monte Carlo data of the CP 3 model. It is 
found that in the non-flattening case, the result of the MEM agrees with that of the Fourier transform, while in 
the flattening case, the MEM gives smooth 2(8). 



1. INTRODUCTION 

It is well known that QCD in principle has a 
9 term. The 9 term is deeply associated with 
non-perturbative properties of QCD at low en- 
ergy and provides us with interesting issues such 
as the strong CP problem, possibilities of rich 
phase structures in 9 space and so on. So it is 
a challenging subject to investigate the dynamics 
of QCD with the 9 term. 

The 9 term in lattice field theory makes the 
Boltzmann weight complex and prevents one from 
performing Monte Carlo (MC) simulations di- 
rectly. This is the sign problem. This problem 
can be circumvented by Fourier-transforming the 
topological charge distribution 
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where Q, (j> and S are the topological charge, a 
field of the system and an action, respectively. 
The measure [2?</>]q implies that the integral is 
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restricted to configurations of <f> with Q. The par- 
tition function is given in terms of P{Q) by 



z{o)=Y J p{Qy 



(2) 



Although this method works well for small vol- 
umes, it does not work for large ones because er- 
rors of P(Q) affect strongly the behavior of the 
free energy density, f(9) = — y\ogZ(9) (V is a 
volume). This is called flattening |4l5| . Flatten- 
ing can be remedied by reducing the errors, but 
this is hopeless, because exponentially increasing 
statistics are needed as volume increases. 

In order to deal with flattening, we have uti- 
lized the maximum entropy method (MEM) |fil7l 
|H]. In ref.JH], we applied the MEM to mock data 
of the Gaussian P(Q) to study whether the MEM 
is effective to our issue. In the non-flattening 
case, the MEM reproduces the exact Z(9), while 
it gives smooth Z(9) in the case where flatten- 
ing occurs in the Fourier method. In the present 
work, as the next step, we apply the MEM to 
MC data of the CP^ -1 model, which has several 
dynamical properties in common with QCD. 
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2. FLATTENING AND MEM 

We simulated the CP 3 model with a fixed point 
action at a fixed coupling constant, /3 = 3.0[3] ■ 
The lattice extension L is changed from 4 to 96. 
The statistics are more than 1 million for each 
case. 
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Figure 1. Free energy density f(9) of the MC 
data of the CP 3 model. Lattice extension L of 
the system is changed from 6 to 56. 



Figure displays f(9) obtained numerically by 
Fourier-transforming P(Q) of the MC data for 
various volumes. Up to L — 38, the Fourier trans- 
form works, but f(9) for L = 50 and 56 cannot be 
calculated properly. Especially for L = 56, f(9) 
becomes flat for 9 > 2.0. This is nothing but flat- 
tening. The density f(9) for L = 50 breaks down 
for 9 > 2.3 due to the negative values of Z(9). 
We also call it flattening because errors of P(Q) 
disturb the behavior of f(9). 

In order to deal with flattening, we employ the 
MEM. To this end, we utilize the inverse Fourier 
transform. The MEM is based on Bayes' theorem 
and gives the most probable image of Z{9) . In our 
case, the most important object is the posterior 
probability piob(2(9)\P(Q), I) which is the prob- 
ability that 2(9) is realized when the MC data of 
{P{Q)} and information / are given. Information 
/ represents our state of knowledge about Z(9). 
Here we impose the criterion Z(9) > 0. 



The probability prob(Z(6»)|P(Q), I) is repre- 
sented in terms of \ 2 an d the entropy S: 



pxoh(2(9)\P(Q),I) oc exp 



(3) 



where a is a real positive parameter. Convention- 
ally the Shannon- Jayncs entropy is employed; 



S = 



(19 



Z{9) - m{9) - Z{9) log 
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where m(9) is called the default model and is cho- 
sen so as to be consistent with information /. To 
sum up, our task is to calculate the most probable 
image 2(9) such that prob(Z|P(Q), /) is maxi- 
mized( see refs.|8!9| for details). 

3. RESULTS 

We apply the MEM to such MC data that flat- 
tening occurs in the Fourier method as well as 
to such MC data that it does not. Here we use 
the MC data for L — 38 (data A) as an example 
of the non-flattening case and for L = 50 (data 
B) for the flattening one. Two types of default 
model are used: (i) Gaussian type, ma(9) = 
exp[— 7 lo ^ 2 10 fl 2 ] , where a parameter 7 is changed 
from 0.6 to 5.0. (ii) weak coupling region type, 
m w (9), which will be explained later. In the anal- 
ysis, the Newton method with quadruple preci- 
sion is used to calculate an image Z(9) with high 
precision. 

Figure[3displays the most probable image Z(9) 
of data A, as an example of the non-flattening 
case. The Gaussian type default models 171g(0) 
with 7 = 0.6 and 1.0 are used. The parti- 
tion function Z(9)(= Zp(9)) obtained by Fourier- 
transforming P(Q) of data A is also displayed. 
The result of the MEM agrees with that of the 
Fourier transform in the whole 9 region. Note 
that errors of 2(6) are too small to be visible, 
where the error of 2(9) means the uncertainty of 
.2(0) 08. 

Next, let us turn to the data B, the flattening 
case. The results are shown in Fig. [3J The par- 
tition function 2p(9) of data B breaks down for 
9 > 2.3 because the values of 2p(9) become neg- 
ative in this region. For the MEM analysis, the 
Gaussian type ma(9) with 7 = 5.0 and the weak 
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Figure 2. The most probable image Z(9) of the 
MC data for L = 38. The Gaussian type default 
models with 7 = 0.6 and 1.0 are used. As a com- 
parison, the result of the Fourier transform(x) is 
also displayed. 



Figure 3. The most probable image Z(9) of the 
MC data for L = 50. Two types of default model 
are used. The result of the Fourier transform( x ) 
is also displayed. 



coupling region type m w (6) are used. Since P{Q) 
of data A was successfully Fourier-transformed to 
Z{9) = Za{0), we use Za{6) as m w (9) for the 
analysis of data for larger volumes. The results 
of the MEM reproduce smooth Z(9) in the whole 
9 region. For 9 > 2.7 they depend on the de- 
fault models and give no definite solution about 
the behavior of Z{9) for 9 ~ ir due to the large 
magnitude of the errors. For 9 < 2.7, however, 
the results of the MEM show no m(#)-dependence 
and their errors are small. We thus obtain a rea- 
sonably good solution of Z{9) for L = 50 up to 
9 = 2.7, which is wider than the valid region of 
the Fourier method. 

4. SUMMARY 

We applied the MEM to the MC data of the 
CP 3 model for various volumes. It was found that 
the MEM is applicable to the CP 3 model in the 
non-flattening and the flattening cases, and that 
the MEM has the advantage of calculating Z(9) 
for somewhat wider 9 region than the Fourier 
transform at least for L = 50. The following sub- 
jects remain to be investigated; (i) a systematic 
study of the m((5)-dependence of Z(9) and its er- 
ror, (ii) investigation on a question whether the 



MEM is more effective than the Fourier method, 
and (iii) a study of the feasibility of distinction 
between flattening and the signal of a first order 
phase transition 10 . 
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